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ABSTRACT 

We report multifrequency observations with the NRAO Very Long Basehne 
Array (VLBA) of the compact radio sources J0128+6306 and J0547+2721, which 
are viewed through the supernova remnants G127. 1+0.5 and S147, respectively. 
Observations were made at frequencies of 1.427, 1.667, 2.271, and 4.987 GHz. 
The hues of sight to these sources pass through the shock wave and upstream 
and downstream turbulent layers of their respective supernova remnants, and 
thus might detect cosmic-ray generated turbulence produced during the Fermi 
acceleration process. For both sources, we detect interstellar scattering, char- 
acterized by a component of the angular size which scales as the square of the 
observing wavelength. The magnitude of the scattering is characterized by an 
effective scattering angular size 9 so at a frequency of 1 GHz of 13.2 ±2.6 milliarc- 
seconds (mas) for J0128-I-6306 and 6.7 ±2.2 mas for J0547±2721. These angular 
sizes are consistent with the "incidental" scattering for any line of sight out of the 
galaxy at similar galactic latitudes and longitudes. There is therefore no evidence 
for enhanced turbulence at these supernova remnants. We establish upper limits 
to the supernova remnant-associated scattering measures of 8.1 — 14.8 m~^°/^-pc 
for J0128±6306 and 3.0 m-20/3_pc for J0547+2721. 

Subject headings: supernova remnants — cosmic rays — plasmas — turbulence — waves 

1. Introduction 

It is widely believed that the cosmic rays are accelerated by supernova remnants (SNR) 
via the Fermi Class I mechanism. The arguments in favor of this viewpoint are (1) the 
Fermi class I mechanism has been demonstrated to be an effective acceleration mechanism 
at shocks in the solar wind (e.g. Lee (1982)), (2) supernova explosions are one of the 
few obvious phenomena which generate sufficient power to account for the cosmic rays, and 
(3) the detection of synchrotron x-rays associated with the supernova remnant SN 1006 
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indicates electron acceleration to extremely high energies (Dyer et al (2001, 2004)). This 

last observation suggests that ions are accelerated to very high energies as well. Nonetheless, 
direct, or at least entirely convincing demonstration of ion acceleration to energies of ~ 
10^^ eV or greater is not available. However, observational evidence for ion acceleration in 
supernova remnants has recently been considerably strengthened by detection of TeV 7-ray 
emission from the SNR G347.3-0.5 (Aharonian et al (2004)). 

The present paper deals with a radioastronomical method of searching for an indicator 
of ion acceleration at supernova remnants. The method is motivated by observed properties 
of shock waves in the interplanetary medium. In the solar wind, ion acceleration at shock 
waves such as the Earth's bow shock or traveling interplanetary shocks is closely connected 
with the generation of large amplitude magnetohydrodynamic waves (Lee (1982)). Exam- 
ples in the case of the Earth's bow shock are shown in Spangler et al (1988), which also 
provides references to the extensive literature. These waves are generated by cosmic rays 
streaming away from the shock, and constitute a "foreshock" upstream of the shock proper. 
These waves serve to isotropize the particles without reducing their energy, so they may be 
repeatedly overtaken by the shock and thus accelerated. 

These waves also generate substantial density fluctuations (Spangler et al (1988); Span- 
gler, Leckband, and Cairns (1997)), which will modify the refractive index for radio waves. 
In the case of the Earth's foreshock, these waves are generated on a scale defined by the 
cyclotron resonance condition for the reflected ions that generate these waves. Since these 
reflected ions have speeds many times the Alfvcn speed, the excited waves have wavelengths 
much greater than the ion inertial length Va/^i, where Va is the Alfven speed and VLi is the 
ion cyclotron frequency. Spangler, Leckband, and Cairns (1997) report that the dominant 
wavelength for the upstream waves is typically 80-160 ion inertial lengths. 

Radio waves propagating through a region of shock-associated plasma waves will be scat- 
tered, producing detectable radio scintillation phenomena, if the region is extensive enough. 
Woo et al (1985) and Woo and Schwenn (1991) reported enhanced radio wave scattering 
associated with interplanetary shock waves. These papers describe numerous cases in which 
strong scattering transient events appeared on lines of sight passing close to the Sun. It is 
possible to attribute these transients to interplanetary shock waves because the transients 
were preceded by a solar flare or coronal mass ejection (CME), and followed by direct de- 
tection of a shock in interplanetary space by a spacecraft such as Helios. It is not clear from 
these observations how much of the enhanced scattering is due to turbulence in the foreshock, 
and how much is due to the downstream shocked plasma or the driver gas associated with 
the CME. All of these regions have counterparts in supernova remnants. 

Radio scintillation observations further out in the interplanetary medium, at heliocentric 
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distances ~ 0.5 a.u., and typically made at longer radio wavelengths, also show scattering 
transients. These latter transients are usually associated with corotating interaction regions 
(CIR) rather than shock waves (Ananthakrishnan, Coles, and Kaufman (1980)). However, 
during particularly large flares or CMEs, transients can be observed at relatively long wave- 
lengths (e.g. 90cm) and at heliocentric distances of several tenths of an a.u. An example 
is the large transient associated with the Bastille Day flare and halo CME in July 2000 
(Tokumaru et al (2003)). Obviously, the most powerful interplanetary shocks provide the 
best analogs for supernova remnant shock waves. 

The idea of looking for phenomena similar to those observed in hcliosphcric shocks at 
SNR shocks was the basis of radioastronomical observations carried out and analysed nearly 
two decades ago by Spangler et al (1986) (Paper 1) and Spangler, Fey, and Cordes (1987) 
(Paper II). In those projects, a very specific prediction of Fermi class I models was examined, 
i.e. that intense turbulence, and thus enhanced radio wave scattering, would be found in the 
foreshock ahead of the main shock jump identified in x-ray or radio continuum data. The 
idea is illustrated in Figure 1, which is taken from Paper II and adapted as described below. 

In Papers I and II, sources were found whose lines of sight passed close to, but upstream 
of the supernova shock as identified on the basis of radio synchrotron radiation. Thus the 
line of sight would not have passed through the downstream shocked region, but would 
have passed through the turbulent foreshock. The idea of these investigations was that 
detection of enhanced radio wave scattering in this foreshock would have been a strong 
indicator of enhanced waves or turbulence in this region, and thus support for the Fermi 
class 1 mechanism. 

No convincing cases were found of enhanced scattering which could be attributed to 
enhanced turbulence in a SNR foreshock. On the contrary, a number of cases were found in 
which a line of sight to an cxtragalactic radio source passed very close to a supernova shock, 
and yet the source showed a very low level of radio wave scattering in the form of angular 
broadening. The best examples were the source 3C418 near the supernova remnant HB21 
(Paper I), the source 0503-1-467 near the remnant HB9 (Papers I, II), and, to a lesser extent, 
the source 0016-1-731 near CTA 1 (Paper I). The most likely explanation for the failure to 
detect enhanced radio wave scattering is that the radio lines of sight, although close on the 
sky to the SNR shells, passed upstream of the turbulent foreshock regions. In terms of the 
parameters defined in Figure 1, this would require the foreshock thickness AR < y where y is 
the closest approach of the line of sight to the SNR shock. The plausibility of this suggestion 
is discussed in Papers I and II. 

In the present paper, a modification of the approach of Papers 1 and 11 is employed, 
which circumvents the possibility that the line of sight could pass upstream of the turbulent 
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foreshock. We now observe radio sources viewed through supernova remnants, with the hne 
of sight indicated by the dotted hne in Figure 1. The advantages of observing a radio source 
through a supernova remnant are twofold. First of all, as is obvious from the cartoon, the line 
of sight must pass through the foreshock, so the radio measurement is affected by this region 
at some level. Second, the mean plasma density is higher in the downstream region by the 
shock compression factor, so there is more plasma to fluctuate. Finally, theoretical arguments 
indicate that shock waves should act as amplifiers of MHD turbulence, in the sense that the 
magnitude of fluctuations in magnetic field and plasma velocity will increase at the shock 
discontinuity (McKenzie and Westphal (1968, 1969); Zhuang and Russell (1982); Vainio 
and Schlickeiser (1998)). Since these magnetic and velocity fluctuations are the drivers of 
the density fluctuations, the density fluctuations responsible for radio wave scattering effects 
should be magnified at the shock as well. Thus, one might expect the SNR shock to increase 
the detectability of the turbulence which is upstream of the shock. 

The obvious disadvantage to downstream observations is that detection of enhanced 
scattering is ambiguous, in that it could be due to turbulence generation mechanisms in 
the post-shock fluid which are unrelated to plasma processes in cosmic ray acceleration. 
Nonetheless, any detection of enhanced radio scattering in the vicinity of a SNR shock 
would be of interest, and the magnitude of the effect can be considered an upper limit to 
turbulence in the foreshock. 

The specific sources which have been chosen for observation arc J0128+6306, an extra- 
galactic radio source which appears nearly dead-center in the supernova remnant G127.1+0.5, 
and J0547+2721, which is viewed through the limb of the radiative supernova remnant S147. 
More description of the sources and associated remnants is given in Section 3 below. 

It might be thought that the enhanced scattering due to a small, highly turbulent 
region in a supernova remnant, characterized by a size scale of a few parsecs, could not 
produce a detectable enhancement of the general interstellar scattering which is due to the 
effect of fines of sight of several kiloparsecs in length. This point is addressed quantitatively 
in Section 5 below. However, an immediate comment is that HIT regions and interstellar 
bubbles associated with star formation regions, which have similar spatial extents, are very 
strong scatterers of radio waves and dominate the "incidental" scattering of the interstellar 
medium (Cordes, Weisberg, and Boriakoff (1985); Spangler and Cordes (1998)). It is 
therefore plausible that highly atypical conditions near a SNR shock, such as high levels of 
wave turbulence, could produce interstellar scattering that dominates that of the rest of the 
interstellar medium. 

The organization of this paper is as follows. Section 2 briefiy summarizes the relevant 
aspects of interstellar scattering for this project. Section 3 describes the details of the 
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observations and the way in which the data were processed to yield the information of 
interest. Section 4 yields the observational results, including our estimates of the degree of 
scattering in the two supernova remnants. Section 5 interprets those results in terms of the 
properties of turbulence near the SNR shocks of our two remnants. Finally, Section 6 gives 
the conclusions. 



2. Relevant Formulcis for Interstellcir Radio Scintillation 

In this paper, we will employ the same techniques and observables as in Papers I and 
II. The purpose of this brief section is to introduce the vocabulary and formulas that will be 
used later, and to clarify the sort of measurements which are made. 

Interstellar turbulence possesses plasma density fluctuations, which subsequently pro- 
duce random radio wave propagation (Rickett (1990)). These density fluctuation are de- 
scribed in terms of their spatial power spectrum, Psn{q), where q is the spatial wavenumber 
of the fluctuation. The standard form assumed for this power spectrum is 

Psniq) = C^g-" (1) 

where Cj^ is the normalization constant of the density power spectrum, and will be the 
primary physical quantity of interest in this paper. In intuitive terms it may be thought of 
as the "intensity" of the turbulence. The spectral index a has been shown to be close to the 
Kolmogorov value of a = 11/3 for many environments in the interstellar medium, and that 
value will be assumed in this investigation. 

Radio waves propagating through the interstellar medium undergo a host of random 
propagation effects, such as intensity scintillations, "corrugation" of the spectrum, and angu- 
lar broadening or blurring. This paper concentrates on measurements of angular broadening, 
which is a phenomenon which can be measured with Very Long Baseline Interferometers. 
We will measure the excess angular size due to scattering in the interstellar medium. 

A point source viewed through a turbulent plasma with a power law spectrum of density 
fluctuations is not precisely a Gaussian image, but the difference between a best-fltting 
Gaussian and the true intensity distribution is relatively subtle and difficult to measure (e.g. 
Spangler and Cordes (1988)). It has been demonstrated that the best- fit Gaussian angular 
size can yield good quantitative measures of turbulent properties (see Cordes, Weisberg, and 
Boriakoff (1985)). This matter is discussed further in Section 5.3 below. 

Radio propagation measurements yield a value for the scattering measure, which is 
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defined as 

SM^ f Clds = ClZ,ff (2) 

J LOS 

where the integral is taken to be along the line of sight. In the second expression, Zg// is 
the effective thickness of a turbulent layer causing the radio propagation effects. 

Extragalactic radio sources, which arc used as the probes in the present investigation as 
well as those of Papers 1 and 11 and many other investigations, are not point sources, and it is 
not always the case that the intrinsic size is much smaller than the angular broadening size. 
In this case, the measured structure of the radio source is the convolution of the intrinsic 
structure with the angular power pattern due to interstellar scattering. When the intrinsic 
structure is barely resolved by the interferometer, and can be modeled as a Gaussian, the 
Gaussian-equivalent angular size of the source as measured is the quadratic sum of the 
intrinsic size {9i) and scattering angular size {Os)- 

The intrinsic and scattered components of the angular size can be separated via multi- 
frequency observations, because they have different dependences on the frequency of observa- 
tion. The scattering size scales as very close to the inverse square of the radio frequency; this 
is a consequence of the exact wavelength-squared dependence of the phase structure function 
of a scattered source. The inverse-frequency-squared dependence of the effective angular size 
has been verified for well-observed and analysed sources such as 2013+370 (Spangler and 
Cordes (1988)). The intrinsic size, on the other hand, has a much less pronounced depen- 
dence on the frequency of observation. For the typical case of a self-absorbed synchrotrion 
source with a flat spectrum, the angular size is approximately inversely proportional to the 
frequency (Marscher (1977)), a theoretical result which is borne out by extensive observa- 
tional experience. 

A model expression for the measured angular size of a source which is affected by both 
interstellar scattering and intrinsic structure is (Fey, Spangler, and Mutel (1989)) 

^o'J") = f + I' (3) 

where Oobsi^) is the measured Gaussian equivalent size at a frequency z/ (in GHz), Ojq is the 
intrinsic angular size at a frequency of 1 GHz, and ^^50 is the scattered angular size at 1 GHz. 
Multifrequency observations can solve or fit for the parameters 9 10 and 6*50- By convention, 
all angular sizes are taken to be full width at half maximum (FWHM) of the brightness 
distribution. 

Once the scattering size 6 so is known, the following formula can be used to obtain the 
scattering measure (equation (2) of Paper I, obtained from Cordes, Weisberg, and Boriakoff 
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(1985)) 



^50 = 2.24 (^ClfZeff) 



3/5 



mas 



(4) 



where Z^ff is in parsecs. 

The analysis employed in this paper uses the concepts and formulas presented above. 
We acknowledge that this is far from the most sophisticated analysis which can be done 
with radio scattering measurements. However, the sort of observations discussed here, which 
were successfully employed in Papers I and II and Fey, Spangler, and Mutel (1989), can be 
considered as analogous to a litmus paper test, i.e. a relatively simple measurement that can 
test for the presence of interstellar scattering, and give a good quantitative measure of its 
magnitude as described by the scattering measure. Values of the scattering measure obtained 
in this way compare very well with subsequent observations utilizing more sophisticated 
analyses of a larger amount of data (Fey, Spangler, and Mutel (1989)). 



Following the negative results of Papers I and II, a literature search was made for 
supernova remnants with strong, compact radio sources in the background, so the line of 
sight to the compact source would have the geometry shown in Figure 1. This search yielded 
two particularly good cases in which a relatively strong, compact source, suitable for VLBI 
observations, was viewed through a supernova remnant. The first is J0128+6306, which 
is seen through the supernova remnant G127.1+0.5, and the other is J0547+2721, viewed 
through the supernova remnant S147. 

J0128+6306/G127.1+0.5 A radio continuum image of the supernova remnant G127.1+0.5 
is given in Joncas, Roger, and Dewdney (1989), who provide a current discussion of this 
SNR. The 1420 MHz continuum image of Joncas et al also clearly shows a bright radio source, 
nearly exactly in the center of this remnant. This is J0128+6306. Although once proposed 
as possible SS433-type object, J0128+6306 has been firmly established as an extragalactic 
source through an optical ID and via HI absorption measurements (Kaplan et al (2004)). 
Kaplan et al (2004) give a distance of 1.3 kpc for G127.1+0.5, based on an association with 
the open star cluster NGC 559. The supernova remnant has apparently not been detected 
in x-rays, so information is not available on the shock speed and density of the upstream 
interstellar medium. 



3. 



Observations and Data Reduction 



3.1. 



Sources Observed 
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J0547+2721/S147 The supernova remnant S147 is a prominent optical SNR. The radio 
structure of this remnant is discussed in Kundu et al (1980). Kaplan et al (2004) give a 
distance of 1.2 kpc, but other estimates range from 800 parsecs to 1.8 kpc. An argument 
in favor of the closer distance is given by the observations of Salmen and Welsh (2004), 
who observed interstellar absorption lines in the spectra of distant stars. They observed 
high velocity absorption lines, obviously associated with the supernova remnant, in stars 
with distances of 880 and 1800 parsecs. The observations of Salmen and Welsh (2004) 
show relatively low ionization state absorption hues out to velocities of ~ ±70 km/sec. This 
observation indicates a relatively slow shock wave with a speed < 100 km/sec, which is 
consistent with other observations cited by Salmen and Welsh (2004), as well as the absence 
of x-ray emission (Sauvageot, Ballet, and Rothenflug (1990)). All the evidence indicates that 
S147 is an old, slow, radiative supernova remnant. The compact radio source J0547+273 is 
located at the extreme edge of S147, just inside a faint outer filament, so the line of sight 
does pass through the interior of the remnant. 



3.2. Observations 

Observations of the radio sources J0128+6306 and J0547+2721 were made with the 
Very Long Baseline Array (VLBA) of the National Radio Astronomy Observatory^ on 6 
October 2002. Observations were also made of the strong compact radio source J0555+3948 
to monitor array performance and guarantee quality of the calibration. 

Observations were made of all three sources at frequencies of 4.987, 2.271, 1.667, and 
1.427 GHz. The intermediate frequency bandwidth was 64 MHz at all frequencies except 
1.427 GHz, where the need to avoid radiofrequency interference limited the recorded band- 
width to 32 MHz. The observational protocol was to observe each of the three sources for 
ten minutes at a given frequency, then tune to the next frequency and repeat a set of ten 
minute scans. The duration of the observing session was slightly over 10 hours. 

The data tapes were correlated with the VLBA correlator in Socorro, New Mexico. The 
correlated data were processed in the AlPS software package. The program CALIB was 
used to fringe-fit the data for fringe rate and delay. Subsequent to removal of the fringe 
rate and delay, the visibility data were averaged over time and frequency to produce one 
complex visibility measurement per baseline per 30 second time interval. The amplitudes 
were calibrated by use of measured system temperatures and known antenna sensitivities. 



^The National Radio Astronomy Observatory is managed by Associated Universities Incorporated under 
contract from the National Science Foundation 
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More limited, but nonetheless useful observations of J0128+6306 and J0547+2721 were 
made with a Mark III VLBI interferometer in 1988, at frequencies of 4.99 and 1.65 GHz. A 
description of the technical details of those observations is given in Section 3.2 of Spangler 
and Cordcs (1998), which used data from the same observing session, and thus is not 
repeated here. 

3.3. Data Analysis 

The primary measurement was a fit to the data in the (u,v) plane for a Gaussian 
equivalent angular size of the source. The correct interpretation of this measurement requires 
that it refer to a compact, structurally simple source. The ideal situation is one in which the 
intrinsic size is unmcasurably small, and all structure is due to interstellar scattering. To 
insure that the analysis was undertaken on isolated, structurally-simple sources, we mapped 
and cleaned all three sources. The maps of J0547-I-2721 showed simple, barely- resolved 
structure at all frequencies, so we could proceed directly with the analysis described in 
Section 2. 

The source JO 128+6306 showed evidence of extended structure at the lower two fre- 
quencies. This is indicated in Figure 2, where we show the maps at (a) 1.427 GHz and (b) 
4.987 GHz. At the higher frequency, the source is simple and barely resolved. At 1.427 
GHz, there is extended emission about the compact component. This is very common with 
such sources, and the extended emission is more conspicuous at lower frequencies due to 
its steeper spectrum. The effect of this emission needs to be removed before an analysis 
employing equations (3) and (4) can be done. 

We therefore used the AIPS task UVSUB to subtract the CLEAN components corre- 
sponding to the extended structure shown in Figure 2 from the 1.427, 1.667 and 2.271 GHz 
visibility data for J0128-I-6306. This corrected (u,v) data set was then fit for an effective 
angular size Ogbs- 

The calibrated (u,v) data (corrected, if need be for extended emission) were fit for 
the effective Gaussian angular size in three ways. First, we wrote a routine in Mathcad to 
produce a least squares fit to the visibility amplitude V{u,v) = Sc{u,v)/ Sq where Sc{u,v) is 
the correlated flux density and 5*0 is the total source flux, to yield a value of 6obs- We also 
used the AIPS task UVFIT, and operated it in two ways. In the flrst, the total flux density 
5*0 (or equivalently, the zero spacing flux) was fixed at its measured value. In the second 
approach, 5*0 was also treated as a parameter to be fit, so it could be adjusted within a small 
range. The full {u, v) range of the measurements was not included in these fits. Baselines 



-10- 



were not included for which the source was significantly resolved (e.g. a visibility less than 
~ 0.20), since for such high degrees of resolution, a Gaussian is not a good approximation 
to the visibility function (see Section 5.3). 

The three methods yielded results in very good agreement. The maximum discrepancy 
of an estimate from the mean of the three estimates is 10 %, which we adopt as the error 
associated with our mcasm^cmcnts. An illustration of one of our (u,v) data sets and the 
best-fit Gaussian visibility function is shown in Figure 3 as the solid line. The dashed and 
dotted lines are discussed in Section 5.3 below. 

The 1988 data were also fit in the (u. v) plane for an effective angular size. In what 
follows, we have only utilized the 1.65 GHz measurements from the 1988 session, since the 
lower frequency gives information on the scattering contribution to the angular size. 

4. Observational Results 

Equivalent Gaussian angular sizes were obtained at all four frequencies for the sources 
J0128+6306, J0547+2721, and J0555+5948. As expected, the effective angular size of the 
calibrator source J0555+5948 was substantially smaller at all frequencies of observation. The 
measured angular sizes for J0128+6306 and J0547+2721 are given in Table 1. The angular 
sizes given in Table 1 were then fit to equation (3), yielding values for the parameters ^50 
and 9 10. Errors on 6*50 and 9 10 were determined by the ranges in each quantity which were 
consistent with an acceptable reduced xl- The hmit which was chosen for the acceptable 
range was xl < 2.60, corresponding to a likehhood of chance occurrence of 5 % for a fit with 
u — 3 degrees of freedom. The values quoted in Table 1 and used throughout the paper are 
confident at the 95% level. 

Values for 9so and 9io for the sources J0128+6306 and J0547+2721 are given in Table 
2. The content of the fourth through sixth columns will be described below. The model 
equation (3) provided a good fit to the data for both sources. This is illustrated in Figure 
4, which shows for each source equation (3) with the angular size parameters given in Table 
2. The fact that equation (3) provides a good fit to the measured angular sizes at the 
four frequencies with the angular size parameters given in Table 2 supports our detection of 
angular broadening due to interstellar turbulence along the line of sight to each of the radio 
sources. 

Interstellar scattering at a measurable level is expected along these low-latitude lines of 
sight, so we next consider how much of the angular broadening illustrated in Figure 4 and 
parameterized in Table 2 is due to the supernova remnants, and how much is "incidental" 
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angular broadening expected for such low- latitude lines of sight, and which would be present 
even in the absence of a supernova remnant. The incidental scattering may be estimated from 
Lazio and Cordcs (1998), who used a large body of pulsar scintillation data to assemble a 
model of galactic radio wave scattering. Lazio and Cordes (1998) reported their estimates of 
angular broadening in terms of our parameter ^^50- The published table in Lazio and Cordes 
(1998) was supplemented by an updated version of the model available on the internet. 
The pubhshed table in Lazio and Cordes (1998) had a value for J0547-I-2721, but not for 
J0128-I-6306. Data are available for both sources in the updated table. The expected values 
of incidental interstellar scattering provided by the Lazio and Cordes model are given by the 
parameter 9inc (column 4 in Table 2). 

Our conclusion from Table 2 is that the measured angular broadening {9so) for our two 
sources is consistent with the expected, incidental broadening for any line of sight through 
this part of the galaxy (^j„c)- Our measured value of ^50 for J0128-I-6306 is marginally larger 
than 9inc, and in fact the nominal value for 6inc is outside (smaller than) the allowed range 
of ^-t the 95 % confidence level. This small excess could justify additional investigation 
of this source- supernova combination, such as by measuring 630 for sources near, but outside 
the SNR G127. 1+0.5. Nonetheless, our value for ^50 only slightly exceeds Oi^c (< 40 %), 
and there is overlap between the statistically-allowed ranges of the two parameters. For the 
remainder of this paper, we will use our measurement of ^50 as an upper limit to scattering 
attributable to the supernova remnant. For the case of J0547-I-2721, the measured interstellar 
scattering is completely consistent with the a-priori estimate of incidental scattering. In fact, 
the nominal value of ^50 is actually less than the nominal value of Bine- 

Another property of radio scattering is anisotropy of the broadened image, a charac- 
teristic which contains information on the anisotropy of the density irregularities and their 
distribution along the line of sight. Image anisotropy, characterized by the axial ratio (AR) 
of the scattered image and the orientation of the major axis of the image, is best measured 
for highly scattered sources whose structure is dominated by interstellar scattering. This 
condition of heavy scattering is not well satisfied for our sources, particularly J0547-I-2721, 
so an extensive investigation of anisotropy is not made here. 

Nonetheless, we did briefiy examine our data for evidence of large image anisotropy 
which might provide qualitative evidence of a localized, anomalous scattering layer. The fits 
to the visibility data described in Section 3.3 returned separate estimates for the major and 
minor axes of the Gaussian brightness distribution (together with errors on these quantities). 
The major and minor axes were combined in the preceding analysis to yield 9obs- The results 
of these fits at 1.427 and 1.667 GHz were examined for evidence of scattering anisotropy. 
At these two lower frequencies the scattering size is maximized relative to the intrinsic size. 
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For J0547+2721 AR = 1.4 ± 0.2 at both 1.427 and 1.667 GHz. Although this value is 
comparable to those reported for sources behind the Cygnus OBI association (Spangler and 
Cordes (1998)), it may also reflect asymmetry of the intrinsic structure. For the more heavily 
scattered source J0128+6306, the axial ratios are 1.20 ± 0.14 at 1.427 GHz and 1.04 ± 0.41 
at 1.67 GHz. 

The main conclusion to be drawn from these numbers is that while anisotropy factors 
between 1.0 and 1.5 might be present for these sources, the data are also consistent with 
isotropic scattering (particularly for J0128+6306). In any case there is no evidence for highly 
anisotropic scattering as might happen if the scattering were dominated by a small region 
with anomalous turbulence along the line of sight. 

In summary, our basic observational result is that there is no significant enhancement 
of the radio wave scattering to these sources due to the supernova remnants. In the next 
section, we provide a quantitative analysis of this result. 



In this section, we use our data to limit the properties of cosmic-ray-generated turbulence 
near supernova remnants. The data in Table 2 can limit the scattering measure associated 
with these two supernova remnants. We make use of formula (4) which relates angular 
broadening to the physical parameter of interest, the scattering measure. 

The scattering measure to a line of sight of interest, containing an object such as a 
supernova remnant which contributes an excess scattering measure, can be written as SM = 
SMi + 6SM. The variable SMi is the incidental scattering measure for that line of sight, 
which would be measured in the absence of the supernova remnant, and 6SM is the excess 
scattering measure due to the supernova remnant. Using equation (4), it is straightforward 
to derive the following formula. 



5. Interpretation of Results 




(5) 



We can calculate SM^ using the estimates of 9inc and equation (4). Since we are reporting 
all our results as upper limits, we likewise will report an upper limit on SSM. 



5.1. J0128+6306 



The most secure upper limit to the interstellar scattering of JO 128+6306 is expressed 
by the upper limit to ^^^o at the 95 % confidence level, 6so < 15.8 mas. Use of this value 
and the value for 9inc given in column 4 of Table 2 gives < 1.34. A value of 9inc = 9.5 
mas corresponds to an incidental scattering measure of 11.1 m^^'^/^-pc. These numbers then 
yield an upper limit to the supernova-remnant-associated scattering measure of SSM < 14.8 
m~^°/^-pc. This value is given in the fifth column of Table 2. 

A smaller, and perhaps more realistic upper limit on SSM (although not the largest 
consistent with the data) is given by using the nominal value of ^50 (i-e. the best-fit value) 
in equation (5), rather than the upper limit at the 95 % confidence level. This value is 13.2 
mas. When this is done, we obtain SSM < 8.1 m~^°/^-pc. This value is listed in parentheses 
in the fifth column of Table 2. 



5.2. J0547+2721 

For this source, the nominal or best-fit value of ^50 actually falls below the estimated 
value of 9inc- This obviously refiects the estimation error associated with both the 9inc and 
^50 values. For purposes of extracting an upper hmit to SSM from these data, we took ^50 
to be the upper limit at the 95 % confidence level. This gave ^50 < 8.9 milliarcseconds, and 
yields a corresponding value of SSM < 3.0 m~^°/^-pc. This value is also given in column 5 
of Table 2. 



5.3. Model Visibility Functions with Meeisured SM 

With these estimates of the scattering measures {SM as well as SSM) for both sources, 
obtained via the methods discussed in Section 2, we can compare our data with a more accu- 
rate model for the visibility functions of these sources. We consider a source which possesses 
intrinsic, resolvable structure that can be approximated with a Gaussian brightness distribu- 
tion. This source is also subject to angular broadening due to Kolmogorov turbulence. The 
visibility function of the source is a product of two functions, one a Gaussian representing 
the intrinsic structure, and the second describing the angular broadening of a point source 
viewed through a medium possessing Kolmogorov turbulence. In the analysis of Section 2, 
this broadening function was also approximated by a Gaussian, but a more accurate function 
is given immediately below. The composite visibility function (the measured correlated fiux 
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density normalized by the total flux of the compact component) is 



V{u,v) = e"4T?2'^?("^+«^)e-5'D^M 



(6) 
(7) 



The function in equation (7) is the phase structure function, and is straightforwardly 
obtained from eq (3) of Spangler and Cordes (1998). In equation (6) 6i represents the 
intrinsic angular size (FWHM) at the frequency of observation, rg is the classical electron 
radius, and f{a) is a constant involving gamma functions of the spectral index a defined in 
equation (1); for a = 3.67, f{a) = 1.18. 

We compared eq (6) with our observations of both sources at 1.427 and 1.667 GHz. At 
these lower frequencies the scattering was most pronounced. The input parameters to eq 
(6) were the intrinsic size 9j, taken from the 9jq values in Table 2 and extrapolated to 1.427 
and 1.667 GHz, and the total scattering measure along the line of sight (SM), which may 
be calculated from 9so in Table 2 and eq (4). For both sources the function (6) provided a 
satisfactory fit to the data at both of the lower two frequencies. This is illustrated in Figure 
3. The sohd curve is the Gaussian model fit to the data for J0128+6306, described in Section 
4 above. The dashed hue shows eq (6) with the best-fit value of the scattering measure. The 
dotted line represents (6) with the upper limit to SM, obtained from the upper limit of 
^50 < 15.8 mas. 

It may be seen that the more accurate visibility function given by eq (6) provides a 
satisfactory representation of the observations and, as expected, better describes the mea- 
surements on long baselines where the source is substantially resolved. This exercise corrob- 
orates the values for the scattering measures SM and SSM obtained with the techniques 
described in Section 2. 



5.4. Comparison with Scattering in the Cygnus OBI Association 

It is instructive to compare the upper limits to the SNR-associated scattering measures 
just quoted with those due to the HII region associated with the Cygnus OBI association 
(Spangler and Cordes (1998)). The most weakly scattered lines of sight there (given in 
Table 2 of Spangler and Cordes (1998)) had scattering measures of at least 140 m~^°/^- 
pc, almost all of which can be interpreted as dSM due to Cygnus OBI structure^. Even 



^The values of scattering measure in Spangler and Cordes (1998) were reported in units of m^^^Z-^-kpc; 
the numbers in Table 2 of Spangler and Cordes (1998) should be multiplied by 1000 to compare with the 
values in the present paper. 
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heavier scattering has been reported for the HIT region NGC 6334A (Trotter, Moran, and 
Rodriguez (1998)). From an observational standpoint, HII regions are more prominent 
scattering entities than supernova remnants. 

5.5. Limits on the turbulence parameter Cj^ 

The local measure of turbulence is the quantity C^^, which is the normalization constant 
of the density spectrum (Eq (1)). This can be obtained from the scattering measurement 
if the thickness of the turbulent region (the parameter Z^jf in equation (4)) is known. As 
discussed in Section 3.1, reasonable distance estimates exist for both supernova remnants, 
being 1300 parsecs for J0128+6306 and 850 parsecs for J0574+2721. The angular radii of 
the remnants can be measured from the radio map in Kundu et al (1980) to be l.°5 for 
S147, and from the map in Joncas, Roger, and Dewdney (1989) for G127. 1+0.5 to be 0°.42. 
We therefore estimate the radii of the supernova remnants to be 9.5 parsecs for G127.1+0.5 
and 22 parsecs for S147. 

For the purposes of this paper, we are interested in the region downstream of the SNR 
shock in which the amplified turbulence from the upstream region remains at high levels. 
This region is analogous to the magnetosheath region of the Earth's bowshock. We do not 
know how extensive this region is. As a rough estimate, we assume that the enhanced 
turbulence fills an annular shell of thickness equal to 10 % of the radius of the SNR shock. 
The total line of sight occupied by this intense "SNR magnetosheath" would then be 20 % 
of the supernova radius. Given the numbers above, Zf.ff would be 4.4 parsecs for S147 and 
1.9 parsecs for G127. 1+0.5. Given the crudeness of our estimate of the size of this region, 
we do not attempt to account for variations in the length of a line of sight through the 
center of the remnant as opposed to one near the limb, as illustrated in Figure 1. In the 
case of G127.1+0.5 this correction would, in any case, be very small since the hue of sight 
goes nearly through the center of the remnant. In the case of S147, the hue of sight to 
J0547+2721 passes near the edge of the remnant, just inside some faint outer filaments. The 
total path length in this case could be larger than twice the radial thickness, so Z^ff could 
be larger, and C% smaller, than we calculate below. However, the precise location of the 
shock front in this part of the remnant is not certain, so we retain the simplest estimate of 
Zf,ff ~ ().20RsNR^ where Rsnr is the radius of the supernova remnant. These calculations 
assume that the bulk of the radio wave scattering occurs in the post-shock, downstream 
region. If turbulence in the foreshock also makes a significant contribution, the true upper 
hmits to C% would be smaller than those we calculate. 

With these admittedly rough estimates of Zg// we can convert the upper limits to 
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scattering measure to upper limits on C^, the mean value of C% in the post-shock region. 
These values are given in the sixth column of Table 2. 

We can again compare these values of for the Cygnus OBI association, given in Table 
2 of Spangler and Cordes (1998). It may be seen that the upper hmits for J0128-I-6306 are 
actually comparable to the measured values for the most weakly scattered Cygnus OBI 
sources, 2013+370 and 2012+383. Although the scattering measures for the Cygnus sources 
are larger, so is the estimated thickness of the turbulent region associated with Cygnus OBI. 
Thus our data cannot exclude the possibility that within the turbulent region of the shock 
in G127. 1+0.5, the turbulent fluctuations in density are as strong as those in parts of the 
Cygnus OBI region. 

Our upper limit for S147 is substantially smaller, < 1 m~^°/^, and smaller than for 
any line of sight through Cygnus OBI. 

Our upper limits to the parameter C^, although small compared to the measured values 
in the ionized shell surrounding the Cygnus OBI association, are very large compared to the 
general background value characteristic of the interstellar medium. The pervasive "Type A" 
turbulence (Cordes, Weisberg, and BoriakofF (1985)), which is associated with the Diffuse 
Ionized Gas (DIG) component of the interstellar medium, has a reasonably uniform value of 

~ 3 X 10-^ m-iV3 (Cordes, Weisberg, and Boriakoff (1985)). Thus, our observations 
are compatible with a factor of between 2000 (S147) to 30000 (G127.1+0.5) enhancement of 
the normal background turbulence level in the turbulent shells of supernova remnants. 

These limits on C]^ can be translated into constraints on the properties of turbulence 
in the foreshock and postshock regions of the supernova remnants G127.1+0.5 and S147, 
as discussed in Papers I and II. Our observations limit a quantity which is determined by 
the rms density fluctuation (T„ and the outer scale to the turbulence Iq (see eq (5) of Paper 
I). Figure 4 of Paper I or Figure 5 of Paper II can be used for an analysis of the present 
observations. This sort of analysis is impeded by our poor knowledge of the plasma density 
either upstream or downstream in these remnants. 

Such an exercise was carried out with our data for J0128+6306 and J0547+2721, using 
the numbers in Table 2 and an equivalent of eq (5) of Paper I. As a representative example 
of such a calculation, we assume that the upstream plasma number density is 1 cm^"^, which 
is increased to 4 cm^^ after passage through the shock. We also adopt, as a representative 
value, an rms density fluctuation which is 18 % of the mean density, so an = 0.72. This 
choice of modulation index is admittedly arbitrary, but is based on observed values in the 
Earth's foreshock (result cited in Spangler (2002)). Our observations then place a lower 
limit on the outer scale of the turbulence of about 0.1 parsecs for J0128+6306. For the 
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more weakly scattered source J0547+2721, the limits are more interesting and restrictive; 
the outer scale would have to be almost as large as the assumed thickness of the turbulent 
shell, Zf.ff/2 — 2.2 parsecs. 

Higher assumed values of Un (resulting from a higher mean density or larger density 
modulation index) make the observational limits more interesting and physically restrictive. 
Such limits could potentially argue against an anomalously elevated level of turbulence in 
the post shock region. On the other hand, a lower value of o"„ than that adopted in this 
calculation is less confining in terms of turbulence characteristics. The situation could be 
improved by x-ray observations of these remnants which would provide knowledge of the 
mean density in the downstream region and the shock speed. 

For the present we are left with the results of this paper, and Papers I and II. There are 
now several cases of radio lines of sight which pass close to, or inside of a supernova remnant 
shock. In no case is there an indication of an enhanced level of plasma turbulence at spatial 
scales of 200 - 10000 kilometers. Since the ion inertial length for a hydrogen plasma with a 
density of 4 cm~^ (our representative value for the downstream density in a SNR shock) is 
~ 120 km, the irregularities to which our measurements are sensitive are from about 2 to 
80 ion inertial lengths in extent. For an upstream region with a density of 1 cm~^ the ion 
inertial length would be larger by a factor of two, and the range of irregularities sampled 
would be approximately 1 to 40 ion inertial lengths. This range of scales is close to, if 
slightly smaller than that on which substantial density fluctuations arc seen in the Earth's 
foreshock. This realization was the motivation for our investigation. However, if physical 
processes in an SNR shock generate compressive irregularities on much larger spatial scales, 
perhaps hundreds to thousands of ion inertial lengths, the anomalous turbulence might have 
escaped detection with our technique. 



6. Conclusions 

The results and conclusions of this paper are as follows. 

1. Observations with the Very Long Basehne Array were made at frequencies of 1.427, 
1.667, 2.271, and 4.987 GHz of the source J0128+6306, which is seen through the 
supernova remnant G127.1+05, and J0547+2721, which is seen through the remnant 
S147. 

2. The measured angular sizes of the compact portions of these sources at the four fre- 
quencies are in good agreement with a model in which the observed source structure 
is a convolution of an interstellar scattering image with the intrinsic structure of the 



-18- 



source. Our estimates of the interstellar scattering size ^50 (FWHM) at 1 GHz is 
13.2 ± 2.6 milliarcseconds (mas) for J0128+6306, and 6.7 ± 2.2 mas for J0547+2721. 

3. The measured values are consistent with a priori estimates of the interstellar scattering 
due to known interstellar turbulence along these lines of sight. We therefore see no 
evidence for enhanced turbulence upstream or downstream of these supernova remnant 
shock waves. Our upper limits on the angular broadening due to SNR-associated turbu- 
lence can be translated to constraints on turbulence at these supernova remnants, but 
the constraints are weak due to poorly known or unknown properties of the supernova 
remnants. 

This work was supported at the University of Iowa by grants ATM99-86887 and ATM03- 
54782 from the National Science Foundation. LGS was also supported by the University of 
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Fig. 1. — Cartoon illustrating the idea of radio probing of the turbulent environment in 
the vicinity of a supernova remnant shock. If the Fermi class I acceleration mechanism is 
operative, there should be a region of intense turbulence upstream of the shock, termed the 
"foreshock". A search for this turbulence may be made by imaging a radio source viewed 
through the foreshock, and see if it is angularly broadened by this turbulence. Investigations 
by Spangler et al (1986) and Spangler, Fey, and Cordes (1987) showed no detectable 
broadening of this sort for several supernova remnants. The present investigation deals 
with a more favorable situation for scattering, in which the line of sight (dotted line) passes 
through the interior of the remnant, where the upstream turbulence is amplified by passage 
through the shock. Figure adapted from Spangler, Fey, and Cordes (1987). 
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Fig. 2.— Maps of the radio source J0128+6306 at (a) 1.427 GHz and (b) 4.987 GHz. The 
lower frequency map shows the extended structure which was removed from the visibihty 
data prior to fitting for an angular size. The beam size is 7.61 by 5.79 milliarcseconds 
at 1.427 GHz and 2.15 by 1.59 mas at 4.987 GHz. The contours in both images are at 
-2, 2, 5, 10, 20, 30, 50, 70, 80, and 90 percent of peak intensity in both images. The peak intensity 
is 0.096 Janskys per beam at 1.427 GHz and 0.168 Janskys per beam at 4.987 GHz. 
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Fig. 3. — Visibility measurements for the source J0128+6306 at a frequency of 1.427 GHz. 
The sohd hne corresponds to the best-fit model visibility function with the Gaussian equiv- 
alent angular size in Table 1. The dashed and dotted lines show a more accurate expression 
for the visibility function of scattered radio source, given by eq (6). The dashed line uses 
the best fit scattering measure SM emergent from our analysis, and the dotted line used the 
upper limit at the 95 % confidence level. 
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Fig. 4. — Observed angular size as function of frequency for (a) JO 128+6306 and (b) 
J0547+2721. The solid curve represents the expected relation for a source whose size is 
a convolution of a scattering angular size with an intrinsic angular size (equation (3)) with 
the angular size parameters given in Table 2. 
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Table 1. Measured Gaussian- Equivalent Angular Sizes 



Source 


Freq (GHz) 


eo6s(FWHM,mas) 


J0128+6306 


4.987 


0.86 ± 0.09 




2.271 


3.26 ± 0.33 




1.667 


4.70 ± 0.47 




1.650 


4.95 ±0.50* 




1.427 


7.23 ± 0.72 


J0547+2721 


4.987 


0.79 ± 0.08 




2.271 


2.13 ± 0.21 




1.667 


2.93 it 0.29 




1.650 


3.40 ± 0.50* 




1.427 


4.69 ± 0.47 



Note. — *: Observations at 1.650 GHz made with 
VLSI network in 1988. 
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Table 2. Retrieved Scattered and Intrinsic Angular Sizes 



Source 


Bio 


Oso 




SSM 


CI 


J0128+6306 


3.4 ± 1.0 


13.2 ±2.6 


9.5 ±1.9 


14.8 (8.1) 


7.8 (4.3) 


J0547+2721 


3.7 ±0.8 


6.7 ±2.2 


9.0 ±1.8 


3.0 


0.68 



Note. — All angular sizes are in milliajrcseconds. The units of the scattering 
measure SM and C% are m ^"/^-parsecs and m 20/3^ respectively. Dual 
values for J0128+6306 correspond to different ways of calculating the upper 
limit to SNR-associated scattering. 



